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ABSTRACT 

The  flow  in  a floating  element  flow  meter  (rotameter)  was  studied 
with  laser  Doppler  anemoraetry  (LDA)  and  computational  fluid 
dynamics  (CFD).  With  the  help  of  refractive  index  matching  to 
eiiminaie  optical  distortion  by  the  flow  tube,  the  laser  Doppler 
anemometer  pemutted  non-intrusive  velocity  measurements  in  all 
regions  of  the  flow  including  the  immediate  vicinity  of  the  float. 
The  computational  investigation  allowed  an  insigh'  to  the  flow 
structure,  particularly  with  respect  to  the  strong  velocity  gradi- 
ents in  the  gap  around  the  float  and  the  recirculation  zone  above 
the  float  Results  were  obtained  for  various  flow  rates,  different 
densities  of  the  float  material  and  two  heights  of  the  float  in  the 
tube 


INTRODUCTION 

Flow  meters  find  widespread  application  for  monitoring  and  centroi 
of  liquid  and  gas  flows  in  mechanical,  autom.otive  and  chemical 
engineering  Many  types  of  flow  meters  rue  in  common  use.  but 
floating  element  flow  meters  have  proven  particularly  appropriate 
for  pipe  diameters  smaller  than  100  mm.  Improvements  ir  such 
meters  necessitate  detailed  investigations  of  the  flow,  using  modern 
techniques  of  experimental  and  computational  fluid  mechanics 

A floating  element  flow  meter  consists  of  a conical  transparent 
vertical  glass  tube  containing  a float  or  •bob’  (figure  1).  In  the 
present  study  the  float  is  formed  of  cylindrical  and  conical  sections, 
although  spherical  floats  may  be  used  in  small  diameter  lubes 
The  float  rises  in  the  tube  until  a balance  is  reached  between  the 
gravitational,  buoyancy  and  drag  forces.  \Vitbin  the  flow  range  uf 
a particular  flow  meter  (depending  on  the  float  shape  and  density, 
the  tube  shape  and  the  fluid  density  and  viscosity)’,  the  height  of 
the  float  in  the  tube  is  linearly  proportioca,  to  the  flow  rate. 

In  early  studies  of  floating  element  flow  meters,  Reppe;  and  Ump- 
fenbach  (1930)  proposed  the  introduction  of  characteristic  dimen- 
sionless quantities,  to  permit  the  use  of  experimentaliy  determined 
drag  coefficients  in  flow  meter  rmalyais  Lutz  (1959)  extended  these 
ideas  by  showing  that  the  transfer  of  flow  coefficients  from  one  flow 
meter  to  another  is  possible  if  geometrical  similarity  exists 

Theories  based  on  similarity  considerations  can  predict  the  varia- 
tion o^fhe  flow  coefficient  with  the  Reynolds  number  in  the  laminar 
and  turbulent  flow  regimes  but  not  in  lam,nar-turbulent  transi- 
tional flow.  Detailed  experimental  and  com  jutational  studies  are 


needed  for  the  transition  region  or  in  regions  where  fast  acceleration 
causes  relaminarization  of  the  flow,  in  order  to  improve  the  accu- 
racy of  floating  element  flow  meters.  For  non-intrusive  velocity 
measurements  inside  the  metering  tube,  a laser  Doppler  anemome- 
ter with  refractive  index  matching  was  used.  A computational  fluid 
dynamics  study  was  made  with  a finite  volume  computer  program 
to  solve  the  Reynolds  equations,  using  a turbulence  model  to  close 
these  equations  so  that  solutions  for  the  entire  mean  flow  field  and 
for  turbulence  properties  could  be  obtained. 

Similarity  analyses  of  the  flow  in  floating  element  flow  meters  are 
reviewed  in  the  next  part  of  this  paper.  The  two  following  sections 
describe  the  experimental  method  and  the  application  of  computa- 
tional fluid  dynamics  to  the  flow.  Measured  and  calculated  velocity 
distributions  are  then  presented  and  discussed 


ANALYSIS  OF  FLOATING 
ELEMENT  PLOW  METER 

The  forces  acting  on  the  float  lead  to  equilibrium  between  the 
weight  of  the  float  pigV^  acting  downwards  luid  the  buoyancy  force 
pg\\  and  the  drag  force  Fj  acting  upwards,  where  Vi  is  the  volume 
and  IS  the  density  of  the  bob,  p is  the  density  of  the  fluid  and  g 
i:  the  gravitational  acceleration: 

PsjV;  = WVi  + F,<  (1) 

The  drag  force  results  from  the  flow  field  surrounding  the  float  and 
depends  e.g  on  the  mean  velocity  profile  of  the  upstrerun  flow  and 
the  separation  in  the  wake  of  the  float.  In  flow  analyses  based  on 
similanty  principles,  these  influences  are  accounted  for  by  empirical 
coefficients  Ci  or  Ct  in  the  drag  law  for 

laminar  flow  rk  = CitiDiU  (2) 

turbulent  flow  Fi  = CrpDlU^  (3) 

where 

o = fluid  viscosity 
Di  = maximum  float  diameter 
U = velocity  in  the  annular  gap  around  the  float  at 
the  mini-mum  cross-section 

The  volume  flow  rate  throug*-  the  rotameter  is 


39.4J. 


Q^^{D^-DhV 

where  D is  Ihe  lube  diarrister  at  the  top  of  the  oob 
Combining  equations  (1)  (2)  and  (4J  gives  for  laminar  flow 


(4), 


(5) 


where  the  constant  q is  defined  as 


rmK 


in  terms  of  the  ot.'.i  ratio  m 
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and  a constant  A'  chtu-acten  tic  of  the  shape  of  the  float 

Under  the  assumption  of  similarity  of  flows  with  the  same  Reynolds 
number,  K can  be  expressed  as 


K = VJDi 

The  Reynolds  number  is  defined  as 

fit  - 


(8), 


(9). 


where  Uin  is  the  velocity  at  the  rotameter  inlet.- 
Using  equation  (3)  instead  of  (2)  y>elds  for  turbulent  flow 


Qt  = 


(10) 


where 


._^rn/K 

^-TV^ 


(11) 


At  low  Reynolds  numbers  the  linear  resistance  law  assumed  in 
equation  (2)  applies  and  q is  a constant  for  given  m.  At  higher 
Reynolds  numbers  the  flow  is  transitional  or  turbulent,  and  logo 
decreases  linearly  with  log  Re. 

Curves  for  0 a'^ainst  Reynolds  number  show  a linear  increase  of 
log  0 with  log  lew  the  laminar  region,  followed  by  a gradual  tran- 
sition to  horizontal  curves  :n  fully  turbulent  flow 


experimental  facility 


Test  Section 


The  experimental  test  section  consisted  of  a glass  flow  meter  *ube 
(Fig  1)  forming  part  of  a closed  circuit  flow  system  with  by-pass, 
through  which  a glycerine  solution  was  circulated  by  a pump  The 
tube  was  mounted  inside  a rectangular  plexiglass  container  with 
two  planar  glass  windows,  which  was  also  filled  with  glycerine.  The 
glycerine  solution  (index  of  refraction  1 455)  minimized  refractive 
index  differences  between  the  fluid  and  the  tube  (index  of  refraction 


1  476),  and  hence  reduced  beam  d.ctortion  due  to  curvaluie  oi  the 
■•ube  wall  The  viscosity  of  the  glycerine  solution  was  strongly 
temperature  dependent,  so  that  the  temperature  of  the  glycerine 
solution  had  to  be  held  consttmt  within  ±0  3 deg  C,  using  a coohng 
circuit  to  remove  the  energy  input  to  the  flow  .-  stem  from  the 
pump 


Fig  1'  Flow  meter  geometry 


Three  different  laminar  flow  cases  were  run  (Table  1),  vtu-ying  the 
flow  rate  (and  hence  Reynolds  number),  the  density  of  the  float 
and  the  height  of  the  float  in  the  tube  (measured  to  the  top  of  the 
float). 


Case  Re  = 4^^ 

1 1S8 

2 400 

3 220 


Float  density 
g/cm^ 
3.85 
7.04 
7 94 


Float  height 
mn. 

160 

160 

105 


Table  1-  Experimental  Cases 


Laser  Doppler  Anemometer 


The  laser  Doppler  anemometer  (Durst  et  al.,  19.‘*1)  was  a forward 
scatter  system  with  a rotating  diffraction  grating  used  for  beam 
splitting  and  frequency  shifting,  similw  to  the  anemometer  desen 
bed  by  Oldengarm  (1976).  The  beam  from  a 15  mW  HeNe  lasei 


was  focused  on  to  the  diffraction  grating  The  diffracted  beams 
of  order  -el  and  —1  were  then  focused  into  the  measuring  volume 
with  the  main  lens,  and  the  zero’th  and  higher  order  diffracted 
beams  were  masked  off  The  half  angle  between  the  intersecting 
beams  .n  the  measuring  volume  was  5.06  deg,  giving  measuring 
volume  dimensions  (1/c^  intensity  level)  of  approximately  0.2  mm 
diameter  and  2 mm  length  The  rotational  speed  of  the  diffraction 
grating  was  adjusted  to  give  a frequency  shift  exceeding  twice  the 
mean  Doppler  frequency 

Forward  scat'ered  light  from  particles  in  the  flowing  glycerine  was 
collected  through  a lens  and  focused  on  to  an  aperture  in  front 
of  a pb'yvomultiplier  The  photomultiplier  signal  was  preamplified, 
band  pass  filtered  and  fed  to  a frequency  tracker  (TSi  model  1096) 
The  analog  output  from  the  tracker  is  a voltage  proportional  to  the 
instantaneous  signal  frequency,  and  hence  to  the  velocity  (with  an 
offset  proportional  t.  the  frequency  shift)  The  output  voltage 
was  low  pass  filtered  and  tent  to  an  ateragng  digtal  voltmeter 
fSolartron  type  JM  1860)  from  which  the  mean  voltage  was  read 

To  minimize  possible  bias  of  the  measurements  resulting  from  drift 
in  the  speed  of  the  rotating  diffraction  grating,  two  measurements 
were  made  at  each  position  with  positive  and  negative  frequency 
sftifts  of  the  same  magnitude  respectively.  The  average  of  the  two 
measurements  (a  voltage  corresponding  to  zero  velocity)  was  then 
subtracted  from  the  voltage  with  positive  frequency  shift,  to  obtrun 
the  velocity  with  the  correct  sign.  The  overall  uncertainty  m mean 
velocity  with  this  procedure  was  estimated  to  be  at  most  2%  This 
accuracy  in  me  laser  Doppler  measurements  is  adequate  to  describe 
the  flow  structure  with  sufficient  detail  for  conclusions  on  improved 
flow  meter  design,  and  to  permit  comparison  with  computational 
results 


COMPUTATIONAL  METHOD 

The  numerical  simulation  of  laminar  flow  m the  floating  element 
flow  meter  vas  carried  out  by  a finite  volume,  multigrid  compu- 
tational program  For  the  mathematical  formulation  of  the  flow 
problem,  the  conservation  equations  for  mass  and  momentum  in  a 
cylindrical  polar  coordinate  system  can  be  written  as 


where  p is  the  density,  ii.t'  and  r,z  are  the  velocity  components 
and  coordinate  directions  as  defined  in  Fig.  1,  p is  the  dynamic 
viscosity  and  P is  the  pressure. 

After  discretization  of  these  partial  differential  equations,  following 
the  method  described  by  Buckle  et  al  (1992),  the  reeulling  coupled 
set  of  equations  for  u,  v and  P is  solved  by  a procedure  based  on  the 
SIMPLE-algoritnm  for  (oiocated  grids  (Penc  et  al , 1988).  After  a 


few  iterations  for  each  variable,  the  coefficients  and  source  terms, 
which  are  functions  of  the  variables,  are  updated  (outer  iterations) 
and  the  whole  process  is  repeated  until  convergence,  as  set  by  a 
prescribed  limit  (typically  10'^  — 10~’)  A “muliigrid  method" 
(Hortmann  et  al ...  1990)  was  used  to  improve  the  rate  of  conver- 
gence of  the  solution.  This  leads  to  a linear  increase  of  computing 
time  with  the  number  of  control  volumes  (CV)  in  comparison  with 
a quadratic  increase  in  the  case  of  single  grid  methods 

The  computational  solution  domain  was  one  half  diametral  plane 
of  the  flow  meter,  and  was  subdivided  into  five  blocks  This  me- 
thod enabled  an  easy  treatment  of  obstacles  inside  the  solution 
domain,  because  the  obstacle  walls  were  identical  with  the  block 
boundaries  Each  block  consisted  of  non-uniform  grids  allowing  a 
better  resolution  near  the  wails. 

The  coarsest  gnd  with  920  CV  Wis  refined  until  a grid  with  14720 
CV  was  reached  Every  CV'  on  the  cc»rge  grid  was  therefore  divided 
into  four  CV  which  form  the  finer  grid. 

The  boundary  conditions  were  u = r = 0 along  the  walls  and 
~ V = 0 along  the  axis  of  symmetry.  At  the  input  boundary  tt-c 
initial  profile  for  the  u-velocity  was  taken  from  the  experiment,  at 
the  outlet  boundary  zero  gradient  was  assumed  for  ail  dependent 
variables  For  case  2 (higher  float  position  at  higher  flow  rate) 
It'  was  necessary  to  extend  the  solution  domain  by  a section  of 
parzdlei-walled  lube  of  diameter  equal  to  the  outlet  diameter  of  the 
flow  meter  lube  and  length  40  mm  in  order  to  obtain  a converged 
solution 


RESULTS 

Measured  Flow  Fielrf 

Velocity  profiles  for  the  three  cases  of  Table  1 were  measured  across 
one  diameter  of  the  tube  at  15  heights.  Fig.  2 shows  results  for 
Case  2 Velocities  are  normalized  by  the  maximum  velocity  at  each 
diameter,  radial  distances  r are  normalized  by  the  tube  radius  R 

The  entry  velocity  profile  is  flatter  than  that  expected  in  fully- 
developed  laminar  pipe  flow,  in  spite  of  the  iov'  Reynold,  number 
The  asymmetry  must  be  attributed  to  the  upstream  piping  At 
the  next  plane  (z  = 50  mm)  the  velocity  profile  is  close’’  to  sym- 
metric and  there  is  a small  reduction  in  the  maximum  velocity. 
Al  z = 90  mm,  however,  a pronounced  asymmetry  appears  again. 
The  peak  velocity  is  now  on  the  other  side  of  the  tube  compared 
with  z = 0.  .suggesting  the  presence  of  some  swirl  in  the  incoming 
flow  Just  upstream  of  ihe  £ (z  = 97  mm)  an  approximately 
symmetric  profile  is  again  found  A depression  on  the  tube  axis 
indicates  an  upstream  propagation  of  flow  disturbance  due  to  the 
float.  In  the  plane  of  the  float  tip  (z  = 100  mm)  deceleration  of 
the  flow  imm.ediately  in  front  of  the  tip  is  very  clear  The  velocity- 
distribution  shows  a cusp  on  the  pipe  axis  with  a local  minimum 
near  0.2  of  the  maximum  velocity.  The  exact  normalized  velocity 
on  the  axis  will  depend  on  how  close  the  measuring  plane  lies  to  the 
plane  of  the  float  tip,  deviations  of  a fraction  of  a mm  will  have  a 
strong  influence  on  the  centre-line  velocity  because  of  a steep  axial 
velocity  gradient  near  the  stagnation  point. 

In  the  next  6 planes  (z  = 120, 134, 143, 151, 152.5  and  160  mm)  the 
progressive  constriction  of  the  flow  in  the  decreasing  gap  between 
the  float  and  the  tube  wall  is  clearly  shown  The  flow  area  is 
reduced  by  s factor  of  2.8  between  the  inlet  station  and  the  plane 
of  the  top  of  the  bob,  while  the  maximum  velocity  increases  only 
by  a factor  of  1.9  between  the  same  stations  The  difference  is 
attributable  to  different  velocity  profiles  in  the  tube  cross-section 
and  in  the  narrow  gap 


A 2 ^ 165  mm  the  flow  undergoes  & rapid  change  as  a result  of 
the  sudden  increase  m cross-sectional  area  immediately  above  the 
float  The  main  flow  separates  from  the  lip  of  the  bob,  leading 
to  a recirculation  zone  behind  the  float.  Between  z = 165  mm 
and  175  mm,  the  peak  re”erse  velocity  increases,  indicating  that 
the  higher  measuring  plane  passes  more  nearly  through  the  centre 
of  the  recirculation  zone  Further  downstream  the  c.^oss-eectiona! 
area  occupied  by  the  main  flow  gradually  increases  and  the  peak 
forward  flow  velocity  falls  correspondingly.  With  increasing  di- 
stance above  the  float,  the  wake  region  indicated  by  a local  ve- 
locity m’.iimu.m  cn  the  tube  axis  steadily  diminishes,  but  at  the 
tube  cutlet  (r  = 210  mm)  there  is  still  a significant  wake  region 
By  contrast,  at  lower  flow  rate  (case  1)  the  effect  of  the  wake  had 
almost  disappeared. 

Comparison  between  experimental 
and  numerical  results 

The  computed  flow  field  for  case  2,  i e.  Reynolds  number  400  and 
the  float  height  160  mm  is  shown  as  a streamline  plot  in  Fig.  3. 
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Fig  3.  Computed  streamlines  (case  2) 
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The  overall  features  of  the  flow  in  the  tube  aie  defined,  and  a 
recirculation  zone  behind  the  “bob"  and  the  reattachment  length 
(i.e.  the  point  where  tt»  — ro-streamhne  intersects  the  symmetry- 
line)  are  shown.  Tue  float  showed  a small  rotation  in  all  flow  cases 
which  was  not  taken  into  account  in  the  numerical  study 
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Fig.  4.  Comparison  of  measure^  rmd  computed  velocity 
profiles  (case  2) 


Experimental  and  computed  velocity  profiles  over  a radius  of  the 
flow  meter  tube  are  shown  in  Fig.  4 for  case  2.  The  comparisons 
a.e  shown  in  selected  measuring  planes  in  the  vicinity  of  the  float, 
where  rapid  streamwise  changes  in  the  velocity  distributions  occur. 
Profiles  are  shown  at  heights  z = 100,  120,  152.5  and  165  mm, 
corresponding  respectively  to  the  tip,  the  conical  portion  and  the 
lower  nm  of  the  bob,  and  a station  just  above  the  bob. 

At  the  first  comparison  station  (z  = 100)  the  computations  indi- 
cate exactly  zero  velocity  on  the  axis,  corresponding  to  a stagnation 
point  at  the  float  tip.  The  measurements  show  a low  but  non-zero 
axial  velocity  at  the  tip.  The  deviation  from  zero  is  probably  at- 
tributable to  two  factors:  slight  unsteadiness  in  the  position  of  the 
float,  even  under  nominally  steady  flow  conditions,  and  smearing 
of  the  radial  distribution  of  velocity  along  the  length  of  the  mea- 
suring volume,  which  was  unavoidably  rather  long  in  companson 
with  the  radius.  This  effect  is  also  apparent  at  radii  up  to  about 
5 mm,  where  the  computations  show  a notably  steeper  velocity 
gradient  than  the  measurements..  The  measured  and  computed 
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peak  velocities  in  the  profiles  agree  rather  well,  rJlhough  the  peak 
of  the  measured  profiV  is  located  at  a lari;er  rfxlius  At  larger  ra- 
dii, where  tlie  velocity  tlcrreases  more  gradually  toward  the  wall, 
computed  ai.d  measured  profiles  are  ir.  close  agreement 

In  the  strongly  converging  annulus  between  the  bob  and  the  tube 
wall,  20  mm  above  the  float  tip  = 120  mm),  the  computations 
reproduce  the  trend  of  the  measureu  res'-lt^  well,  although  there 
are  some  noteworthy  dev  lations  The  measured  points  in  the  region 
of  rising  velocity  adjacent  to  the  bob  are  displaced  to  higher  radius 
relative  to  the  computed  velocity  profiles  In  the  region  of  falling 
velocity  near  the  flow  meter  wall,  however,-  computed  and  mea- 
sured re.sults  deviate  more  noticeably  The  di.screpancy  is  linked 
with  ail  apparent  deviation  from  mass  conservation  in  the  mea- 
surements, which  is  partly  attributable  in  flow  asymmetry  The 
profiles  in  Fig  2 indicate  asymmetry  in  axi«j  velocity  at  equal  ra- 
dii on  either  side  of  the  axis  up  to  10%.  Possible  asymmetry  along 
diameters  inclined  to  the  mam  measurement  planes  r ould  also  have 
contributed  to  the  mass  imbalance 

At  r = 152  5 mm.  measured  and  computed  velocities  apoear  to  be 
in  close  agreement,  except  for  some  deviation  in  the  peak  velocities 
The  confinement  of  the  flow  in  the  narrow  annular  gap  in  this  plane, 
however,  makes  any  differences  difficult  to  resolve. 

5 mm  above  the  float  (z  = 165  mm)  the  flow  is  characterized  bv 
a strong  upwaid  flow  in  an  annular  region  near  the  flow  meter 
tube  wall  and  a recirculation  region  around  the  axis  occupying 
almost  half  the  tube  cross-section  The  forward  velocities  show  a 
reasonably  good  match  between  compulations  and  measurenients, 
but  there  is  a very  marked  discrepancy  in  the  recirculation  zone, 
where  the  measured  nr.aximum  reverse  velocity  is  only  one-third  of 
the  computed  value  and  the  measured  recirculation  region  extends 
radially  to  only  9 mm  compared  with  il  mm  computationally 
Between  9 mm  and  14  mm  in  the  measurements,  however,  is  an 
area  of  almost  zero  mean  velocity,  so  that  the  region  of  forward 
velocity  is  actually  smaller  than  in  the  computations. 

As  a whole,  good  agreement  is  obtained  between  experimental  and 
nume-ical  results,  apart  from  the  separation  region  after  the  float 
Efftcls  caused  by  rotation  of  the  floating  element  and  asymn.ei.y 
of  the  oncoming  velocity  profile  could  not  be  taken  into  .ccounl 
in  t.ie  axi-symmetric  flow  p edictions  A possible  further  enuse  of 
discrepancies  is  an  insuiTicient  number  o!  numerical  nod»-s  down- 
stream of  the  float,  where  the  shear  layer  develops  from  the  upper 
edge  of  t.be  floating  element 

CONCLUSIONS 

This  paper  emerges  from  a sequence  of  studies  at  the  authors’  in- 
stitute related  to  impruvemet-ts  of  voiv.ne  tuid  mass  flow  meters  by 
modern  experimental,  analytics!  and  numerical  techniques  of  fluid 
mechanics.  The  complementaiv  ,tse  of  these  techniques  provides 
the  best  insight  into  the  detail  of  tae  flow  and  their  influence  on 
the  accuracy  of  volume  or  trass  flow  metcr.ng. 

The  results  show  that  the  computational  fluid  dynamics  (CFD) 
code  IS  capable  of  reproducing  the  features  of  the  flow  indicated  by 
the  measurements  by  laser  Doppler  anemometry  Quantitative  di- 
sagreement between  results  by  the  two  methods  does  not  invalidate 
this  conclusion,  but  points  up  the  difficulties  of  Ifiu  "x^ienmental 
and  computational  tasks,  particularly  in  the  narrow  ji  r.o'ar  gap 
existing  in  some  geometrical  planes  T te  result?  give  a eo'-vincmg 
demonstration  of  the  utility  of  modern  tools  of  cornpuiai.unaJ  and 
experimental  fluid  meebames  to  improving  the  analyai.s  and  design 
of  flow  meters  which  traditionally  has  had  to  be  treateo  stm.- 
empirically.  A useful  extension  of  the  work  would  be  a s-udy  for 
Reynolds  numbers  appropriate  to  turbulent  flow  m the  flow  meter 
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